Scaling Dynamics of a Massive Piston in a Cube 
Filled With Ideal Gas: Exact Results 

N. Chernov^''^, J. L. Lebowitz^''*, and Ya. Sinai^ 
February 1, 2008 



Dedicated to Robert Dorfman on the occasion of his 65^^ birthday 

Abstract 

We continue the study of the time evolution of a system consisting of a piston 
in a cubical container of large size L filled with an ideal gas. The piston has 
mass M ^ L"^ and undergoes elastic collisions with ~ gas particles of mass 
m. In a previous paper, Lebowitz, Piasecki and Sinai |LPS| considered a scaling 



regime, with time and space scaled by L, in which they argued heuristically that the 
motion of the piston and the one particle distribution of the gas satisfy autonomous 
coupled differential equations. Here we state exact results and sketch proofs for 
this behavior. 



1 The model and main results 



This paper is a continuation of |P^PS|| , where deterministic scaled equations describing the 
dynamics of a massive piston in a cubical container filled with ideal gas were given. Here 
we state exact conditions on the validity of those equations and outline the arguments. 



Full proofs will be published in a separate paper ||CLl5| ]. 



We refer the reader to [[LP^j and to [pLSj as well as 0, |GPl [KBMj for a 



detailed description of the problem of a massive piston moving in a cylinder. Here we 
just recall necessary facts. 

Consider a cubical domain of size L separated into two parts by a wall (piston), 
which can move freely without friction inside Al. Each part of A/, is filled by a noninter- 
acting gas of particles, each of mass m. The piston has mass M = Mi and moves along 
the X-axis under the action of elastic collisions with the particles. The size L of the cube 
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is a large parameter of the model, and we are interested in the limit behavior as L — oo. 
The mass m of gas particles is fixed. We will assume that M grows proportionally to 
and the number of gas particles is proportional to L^, while their velocities remain of 
order one. 

The position of the piston at time t is specified by a single coordinate X = 
< X < L, its velocity is then given by V" = VL(t) = XL(t). Since the components of 
the particle velocities perpendicular to the x-axis play no role in the dynamics, we may 
assume that each particle has only one coordinate, x, and one component of velocity, v, 
directed along the x-axis. 

When a particle with velocity v hits the piston with velocity V, their velocities after 
the collision, v' and V, respectively, are given by 

V' = {l-e)V + ev (1.1) 

v' = -{1 - e)v + {2 - e)V (1.2) 

where e = 2m/ (M + m). We assume that M + m = 2mL^/a, where a > is a constant, 
so that 

2m a 

When a particle collides with a wall at x = or x = L, its velocity just changes sign. 

The evolution of the system is completely deterministic, but one needs to specify 
initial conditions. We shall assume that the piston starts at the midpoint Xl{0) = L/2 
with zero velocity Vl{0) = 0. The initial configuration of gas particles is chosen at random 
as a realization of a (two-dimensional) Poisson process on the (x, f )-plane (restricted to 
< X < L) with density L^pl{x,v), where pl{x,v) is a function satisfying certain 
conditions, see below, and the factor of is the cross-sectional area of the container. In 
other words, for any domain D C [0,L] x the number of gas particles (x,f) G -D at 
time t = has a Poisson distribution with parameter = L"^ J Jj^ Pl{x, v) dx dv. 

Let Ql denote the space of all possible configurations of gas particles in A^. For each 
realization u E Ql the piston trajectory will be denoted by Xi(t, u;) and its velocity by 

As L — >■ oo, space and time are rescaled as 

y = x/L and r = t/L. (1.4) 
which is a typical rescaling for the hydrodynamic limit transition (see [[LPSj , PLS|| for 



motivation and physical discussion). We call y and r the macroscopic ("slow") variable, 
as opposed to the original microscopic ( "fast" ) x and t. Now let 

Yl{t,u) = Xl{tL,u)/L, Wl{t,u) = Vl{tL,uj) (1.5) 

be the position and velocity of the piston in the macroscopic variables. The initial density 
Pl{x, v) satisfies 

Pl{x,v) = Tro{x/L,v) 
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where the function Tro{y,v) is independent of L. Without loss of generahty, assume that 
ttq is normahzed so that 

r-l roo 

7ro{y,v) dvdy = 1 





Then the mean number of particles in the entire container A/, is equal exactly to 



rL roo 

E{N) = / L^pL{x,v)dvdx = 

Jo J-oo 



In order to describe the dynamics by differential equations, we assume that the func- 
tion 7ro(?/,f) satisfies several technical requirements stated below. 

(PI) Smoothness. no{y,v) is a piecewise function with uniformly bounded partial 
derivatives, i.e. \dno/dy\ < Di and Idno/dvl < Di for some Di > 0. 

(P2) Discontinuity lines. 7ro(?/, f) may be discontinuous on the line y = 1l(0) (i.e., "on 
the piston"). In addition, it may have a finite number (< Ki) of other discontinuity 
lines in the {y, t')-plane with strictly positive slopes (each line is given by an equation 
V = f{y) where f{y) is and < ci < f'{y) < C2 < oo). 

(P3) Density bounds. Let 

7ro(?/,t;) > TTmin > for Vi<\v\<V2 (1.6) 

for some < fi < < oo, and 

sup 7lo{y, V) = VTmax < OO (1.7) 

The requirements (|1.6|) and ( |1.7| ) basically mean that TTo{y,v) takes values of order 
one. 

(P4) Velocity "cutoff". Let 

no{y,v)=0, if |f | < Wmin or |w| > fmax (1-8) 

with some < fmin < v^^ax < oo. This means that the speed of gas particles is 
bounded from above by Vmax and from below by v^^^^. 

(P5) Approximate pressure balance. vro(?/, v) must be nearly symmetric about the piston, 
i.e. 

ko(z/,^^) - vro(l - -t;)| < £0 (1-9) 
for all < ?/ < 1 and some sufficiently small £0 > 0. 
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The requirements (P4) and (P5) are made to ensure that the piston velocity \VL(t, u!)\ 
will be smaller than the minimum speed of the particles, with probability close to one, 
for times t = 0{L). Such assumptions were first made in ||LPS |. 



We think of Di, Ki, ci, C2, vi, V2, fmin, fmax, TTmin and TTmax in (P1)-(P4) as fixed 
(global) constants and Eq in (P5) as an adjustable small parameter. We will assume 
throughout the paper that Eq is small enough, meaning that 

, Ci, C2, Vi, V2y '^min; ^max) ^miii) ^maxj 

It is important to note that the hydrodynamic limit does not require that eq ^ 0. The 
parameter Eq stays positive and fixed as L ^ 00. 
Here is our main result: 

Theorem 1.1 There is an L-independent function Y{t) defined for all r > and a 
positive ~ 2/t>max (actually, r* — > 2/t>max as Eq ^ 0), such that 

sup \Yl{t,uj)-Y{t)\^0 (1.10) 

0<T<r* 



and 



sup \Wl{t,uj) -W{t)\^0 (1.11) 

0<T<T, 



in probability, as L ^ 00. Here W{t) = Y{t). 

This theorem establishes the convergence in probability of the random functions 
Yl{t,uj), W{t,u!) characterizing the mechanical evolution of the piston to the determin- 
istic functions Y{t), W{t), in the hydrodynamic limit L ^ 00. 

The functions Y{t) and W{t) satisfy certain (Euler-type) differential equations stated 
in the next section. Those equations have solutions for all r > 0, but we can only 
guarantee the convergence ( p..lO|) and (|1.11| ) for r < r*. What happens for r > r, 



especially as r — > cxo, remains an open problem. Some experimental results and heuristic 
observations in this direction are presented in ||CL|| and discussed in Section 4. 

Remark 1. The function Y{t) is at least and, furthermore, piecewise on the 
interval (0, r*). Its first derivative W = Y (velocity) and its second derivative A = Y 
(acceleration) remain Eo-smglh sup^|H^(r)| <const-eo and sup^|A(r)| <const-£:o, see 
the next section. 



Remark 2. We also estimate the speed of convergence in (|1.10|) and (|1.11| ): there is 



a Ti > (ri ^ V-Wmax) such that \Yl{t,uj) - Y{t)\ = 0{\nL/L) for < r < n and 
\Yl(t,uj) — Y{t)\ = 0(lnL/L^/'') for ri < r < r*. The same bounds are valid for 
\Wl{t,u!) — W{t)\. These estimates hold with "overwhelming" probability, specifically 
they hold for all G fi^ C such that P{ni) = 1 - 0(L-i°^). 
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2 Hydro dynamical equations 

The equations describing the deterministic function Y{t) involve another deterministic 
function - the density of the gas 7i{y,v,T). Initially, vr(|/,f, 0) = 7To{y,v), and for r > 
the density 7r(y,f,r) evolves according to the following rules. 

(HI) Free motion. Inside the container the density satisfies the standard continuity 
equation for a noninteracting particle system without external forces: 

(^ + ^^) ^(^'^''^) = (2.1) 
for all y except y = 0, y = 1 and y = Y{t). 
Equation ( |2.1| ) has a simple solution 

TT{y,V,T) = TT{y -VS,V,T - s) (2.2) 



for < s < r such that y — vr ^ {0, Y{t — r), 1} for all r G (0, s). Equation has one 
advantage over (^.11): it applies to all points {y,v), including those where the function tt 
is not differentiable. 

(112) Collisions with the walls. At the walls y = and y = 1 we have 

7r(0,w,r) = 7r(0,-w,r) (2.3) 
7ril,v,T) = nil,-v,r) (2.4) 

(113) Collisions with the piston. At the piston y = Y{t) we have 

7r(r (r) - 0, i;, r) = n(Y{T) - 0,2W{t) - v,t) for v<W{t) 
n{Y{T) + 0,v,t) = 7r{Y{r) + 0,2W{t) -v,t) for v>W{t) (2.5) 

where v represents the velocity after the collision and 2W{t) — v that before the 
collision; here 

W{t) = ^Y{t) (2.6) 
dr 

is the (deterministic) velocity of the piston. 

It remains to describe the evolution of W{t). Suppose the piston's position at time 
r is y and its velocity W. The piston is affected by the particles {y, v) hitting it from 
the right (such that y = Y + and v < W) and from the left (such that y = Y — 
and V > W). Accordingly, we define the density of the particles colliding with the piston 
( "density on the piston" ) by 

,(v,r;Y,W) = \^Z^'^'^'^\ < Z (2-7) 

' ' ' ^ 1 7r(F-0,f,r) if V > W ^ ^ 
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(H4) Piston's velocity. The velocity W = W{t) of the piston must satisfy the equation 

/oo 
{v - Wf sgn{v - W) q{v, r; Y, W) dv = (2.8) 
-oo 



We also remark that for r > 0, when ( p.5|) holds, 

/ VTr(Y — 0, V, t) dv J vn^Y + 0, v, r) dv 



W{t) 



j7i{Y -0,v,T)dv j7r{Y + 0,v,T)dv 



i.e. the piston's velocity is the average of the nearby particle velocities on each side. 

The system of (hydrodynamical) equations (H1)-(H4) is now closed and, given ap- 
propriate initial conditions, should completely determine the functions Y{t), W{t) and 
p{y,v,T) for r > 0. To specify the initial conditions, we set p{y,v,0) = Tr{y,v) and 
Y{0) = 0.5. The initial velocity W{0) does not have to be specified, it comes "for free" 
as the solution of the equation (|2.8|) at time r = 0. It is easy to check that the initial 
speed |W^(0)| will be smaller than fmin, in fact W{0) ^ as ^ in (P5). Note that 
if the initial conditions at r = do not satisfy ( p.3|) -( pl5D , there will be a discontinuity 
in p as r — >■ (see also Remark 4 below). 

Equation (|2.8| ) has a unique solution W as long as the piston interacts with some gas 
particles on both sides, i.e. as long as 

mi{v : n{Y + 0, v, r) > 0} < sup{v : n{Y - 0, v, r) > 0} 



Indeed, the left hand side of ( p.8| ) is a continuous and strictly monotonically decreasing 
function of W, and it takes both positive and negative values. The solution W{t) may 
not be continuous in r, though. But if n^y, v, r) is piecewise and has a finite number 
of discontinuity lines with positive slopes (as we require of no{y,v) in Section p, then 
W{t) will be continuous and piecewise differentiable. 

Remark 3. One can easily check that the total mass Ai = JJ TT{y,v,T) dv dy and 
the total kinetic energy 2E = JJ v'^'n'{y,v,T) dv dy remain constant along any solution 
of our system of equations (H1)-(H4). Also, the mass in the left and right part of A/, 
separately remains constant. Equation ( p.8[ ) also preserves the total momentum of the 
gas / / v7i{y,v,T) dvdy, but it changes due to collisions with the walls. 

Remark 4. Previously, Lebowitz, Piasecki and Sinai |[LPS|| studied the piston dy- 
namics under essentially the same initial conditions as our (P1)-(P5). They argued 
heuristically that the piston dynamics could be approximated by certain deterministic 
equations in the original (microscopic) variables x and t. The deterministic equations 
found in [[LPS|] correspond to our (p^ ) -(p^ ) with obvious transformation back to the 
variables x,t, but our main equation (|2.8|) has a different counterpart in the context of 
| LPS|| , which reads 



oo pV 

{v - V{t)fn{Y - 0, V, t)dv- {v- V{t)Y7i(Y + 0, v, t) dv 

V J— oo 



(2.9) 
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Here X = X(t) and V = V(t) = X(t) denote the deterministic position and velocity of 
the piston and n{x,v,t) the density of the gas (the constant a appeared in (|1.3| )). We 
refer to | LPS | for more details and a heuristic derivation of ( p.9| ). Since ( |2.9|) , unlike our 
( pTsD , is a differential equation, the initial velocity V^(0) has to be specified separately, 
and it is customary to set ^(0) = 0. Alternatively, one can set ^(0) = Vl^(O), see ||CL|j| ]. 



Equation ( |2.9| ) can be reduced to (|2.8| ) in the limit L ^ oo as follows. One can show (we 
omit details) that (|2.9| ) is a dissipative equation whose solution with any (small enough) 
initial condition V{0) converges to the solution of ( |2.8| ) during a t-time interval of length 
~ InL. That interval has length ~ InL on the r axis, and so it vanishes as L — >■ oo, 
this is why we replace (p.9|) with ( |2.8|) and ignore the initial condition V^(0) when working 
with the thermodynamic variables r and y. The equation ( p.9|) is not used in this paper. 



We now describe the solution of our equations (H1)-(H4) in more detail. Assume that 
for some r > the gas density 7r(?/,f,r) satisfies the same requirements (P1)-(P4) as 
those imposed on the initial function 7ro(y,f) in Sect. |l|, with constants D[, K[, c[, c'2, 
v[, v'2, fmin5 "^max' ^min vT^ax' whose valucs are uot csseutial, but are independent of 
r. 

We also assume an analogue of (P5), but this one is not so straightforward, since the 
piston does not have to stay at the middle point y = 0.5 at any time r > 0. We require 
that 

|F(r)-0.5|<4 (2.10) 

and for any point {y,v) with v'^^^ < \v\ < v'^^^ there is another point {y^,v^) "across the 
piston" (i.e. such that {y — F(r))(?/* — Y{t)) < 0) satisfying 

ly + y* - 1| < 4, \v + v^\<e'Q (2.11) 

and 

k(l/,t;,r) - 7r(t/*,f*,r)| < 4 (2.12) 

for some sufficiently small e'q > 0. Actually, the map {y,v) \—>- (?/^,,i;^,), which we denote 
by Rr, is one-to-one and will be explicitly constructed below. The constant e'q here, just 
like ^0 ill (P5), is assumed to be small enough, and moreover 

e'o < C^eo (2.13) 

with some constant Cq > 0. 

We now derive elementary but important consequences of the above assumptions. 
Since the density 7r(?/, t>,r) vanishes for \v\ < fmin; so does the function q{v,T;Y,W) 
defined by (|2.7|) . Moreover, for all \W\ < v'^^^, the function q{v, r; Y, W) will be indepen- 
dent of W, and so we can write it as q{v,T;Y). Also, equation ( p^.S] ) can be simplified: 
the factor sgn(f — W) can be replaced by sgnv. Then, expanding the square in (|2.8| ) 
reduces it to a quadratic equation for W: 

QoW^ - 2QiW + Q2 = (2.14) 
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where 

Qo = J sgnv ■ q{v,T;Y) dv (2.15) 

Qi = J V sgnv ■ q{v,T;Y) dv (2.16) 

Q2 = j v"^ sgnv ■ q{v,T;Y)dv (2.17) 
with Y = Y{t). The integrals Qo, Qi, Q2 have the following physical meaning: 

rriQo = mL-mR 

mQ2 = 2{eL - cr) 

where mL,PL,GL represent the total mass, momentum and energy of the incoming gas 
particles (per unit length) on the left hand side of the piston, and m/j,p^, e/j - those on 
the right hand side of it. The value Q2 also represents the net pressure exerted on the 
piston by the gas as if the piston did not move. Of course, if Q2{t) = 0, then we must 



have W{t) =0, which agrees with ( |2.14|) 



Next, under the above requirements on 7r(?/, v, r), the function g(f , r; Y) is, in a certain 
sense, nearly symmetric in v about 1; = (see [PLS | for details). This fact implies that 
Qq and Q2 are small, more precisely 

max{|Qo|,|Q2|}<C"£o (2.18) 

where C" > is a constant depending on the parameters D'^,K[, etc., but not on Eq. At 
the same time, the assumption (P3) guarantees that 

Ql > Ql.min > (2.19) 

where Qi,min is a constant depending on Tr^in, t"^, fg, etc., but not on Eq. 

If Eq is small enough, there is a unique root of the quadratic polynomial ( |2.14| ) on the 
interval (— fmin, "y^in); which corresponds to the only solution of (^]8|). Since this root is 
smaller, in absolute value, than the other root of ( 2.14|) , it can be expressed by 



, X Ql ~ vQi ~ Q0Q2 
W{t) = ^ — (2.20) 

Qo 

where the sign before the radical is "— ", not "+". Of course, ( p.20| ) applies whenever 
Qo 7^ 0, while for Qo = we simply have W{t) = Q2/'2Qi. 

Eqs. ( [^.l^ )-( P^ ) imply an upper bound on the piston velocity: |I^(r)| < B'eq for 
some constant B' > depending on D[, K[, etc., but not on Eq. A similar bound holds 
for the piston acceleration A{t) = W{t), since 

. , , {dQo/dT)W^ - 2{dQ^/dT)W + {dQ2/dT) 
A[t) = 



2{Qi-QoW) 
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and \dQildT\ = \{dQi/dY)W\ <const-eo, see [ pLS|| for more details. 

Next we consider the evolution of a point {y, v) in the domain G := {{y,v) : < y < 
1} under the rules (H1)-(H3), i.e. as it moves freely with constant velocity and collides 
elastically with the walls and the piston. Denote by {yr,Vr) its position and velocity at 
time r > 0. Then (HI) translates into i/r = Vr and Vr = whenever yr ^ {0, 1,Y{t)}, 
(H2) becomes (?/t-+o, fr+o) = {Vr-o, —Vt-o) whenever y^-^Q G {0, 1}, and (H3) gives 

{yr+O, Vr+o) = (l/r-0, 2W (t) - Vr^o) (2.21) 

whenever yr-o = ^(t). Note that ( |2.21| ) corresponds to a special case of the mechanical 
collision rules (|1.1| )- (|1.2| ) with e = (equivalently, m = 0). Hence the point {y,v) moves 
in G as if it was a gas particle with zero mass. 

The motion of points in {y, v) is described by a one-parameter family of transfor- 
mations F'^ : G ^ G defined by F'^{yQ,Vo) = (y^,Vr) for r > 0. We will also write 
F~^{.yT,Vr) = {yo,Vo). According to (H1)-(H3), the density TT{y,v,T) satisfies a simple 
equation 

TT{yr, Vr, t) = n{F''' (yr, Vr), 0) = TToiyo, Vq) (2.22) 

for all r > 0. Also, it is easy to see that for each r > the map F'^ is one-to-one and 
preserves area, i.e. det \DF'^{y,v)\ = 1. 

Now, because of (P4), the initial density tiq{ii,v) can only be positive in the region 

G^ := {{y,v) : < y < 1, Vmin < \v\ < l^max} 

hence we will restrict ourselves to points {y, v) G G^ only. At any time r > 0, the 
images of those points will be confined to the region G'^{t) := F'^{G~^). In particular, 
7r(l/, v,t) = for {y, v) ^ G+(r). 

The map Rr : {y,v) ^— (?/*,f*) involved in (|2.11|) and (|2.12|) can now be defined as 
Rr = F'^ o Rqo F~'^ , where Ro{y,v) = {l—y, —v) is a simple reflection "across the piston" 
at time r = 0. 

We now make an important observation. If a fast point {yr, Vr) collides with a slow 
piston, |W^(r)| <^ \vr\, they cannot recoUide too soon: the point must travel to a wall, 
bounce off it, and then travel back to the piston before it hits it again. 

Therefore, as long as (P1)-(P4) hold, the collisions of each moving point [yr, Vr) G 
G~^{t) with the piston occur at well separated time moments, which allows us to effec- 
tively count them. For {x, v) G G"*" 

N{y, V, t) = #{s G (0, r) : y, = F (s), v^ ^ W{s)} 

is the number of collisions of the point {y,v) with the piston during the interval (0,r). 
For each r > 0, we partition the region G^{t) into subregions 

GUr):={F-{y,v): {y,v)eG+ & N{y,v,T) = n} 

so G^{t) is occupied by the points that at time r have experienced exactly n collisions 
with the piston during the interval (0, r). 
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Now, for each n > 1 we define r„ > to be the first time when a point {yr-, Vr) G G^{t) 
experiences its (n + l)-st collision with the piston, i.e. 

= sup{r > : G++i(r) = 0} 

In particular, ti > is the earliest time when a point (y^jV^-) G G^{t) experiences its 
first recoUision with the piston. Hence, no recollisions occur on the interval [0,ri), and 
we call it the zero-recollision interval. Similarly, on the interval (ri,r2) no more than 
one recoUision with the piston is possible for any point, and we call it the one-recollision 
interval. 

The time moment r^, mentioned in Theorem |TTl| is the earliest time when a point 
{y^^Vr) € G^{t) either experiences its third collision with the piston or has its second 
collision with the piston given that the first one occurred after ri. Hence, r* < r2, and 
actually r* is very close to r2, see below. 

The following theorem summarizes the properties of the solutions of the hydrody- 
namical equations (H1)-(H4). 

Theorem 2.1 ( [|CLS |) Let T > be given. If the initial density TTQ{y,v) satisfies (Pl)- 
(P5) with a sufficiently small sq, then 

(a) the solution of our hydrodynamical equations (H1)-(H4) exists and is unique on the 
interval (0, T); 

(b) the density 7r(?/, f,r) satisfies conditions similar to (P1)-(P4) for all < t < T , it 
also satisfies ( ^.1(\ )- ^J^ ); 

(c) The piston velocity and acceleration remain small, \W{t)\ = 0{eo) and \A{t)\ = 
0{eo); 

(d) we have \Tk — k/v^s,x\ = 0{eo) for all 1 < k < Tvmax, o-nd if Tv^aa^ > 2, then also 
\n - = 0{£o). 



Corollary 2.2 If Eq = 0, so that the initial density TTo{y,v) is completely symmetric 
about the piston, the solution is trivial: Y{t) = 0.5 and W{t) = for all r > 0. 

Lastly, we demonstrate the reason for our assumption that all the discontinuity curves 
of the initial density Tioi^y, v) must have positive slopes. It would be quite tempting to let 
no{y, v) have more general discontinuity lines, e.g. allow it be smooth for fmin < \v\ < fmax 
and abruptly drop to at v = fmin and v = fmax- The following example shows why this 
is not acceptable. 

Example. Suppose the initial density TTo{y,v) has a horizontal discontinuity line v = Vq 
(say, Vq = fjnin or Vq = fmax)- After ouc interaction with the piston the image of this 



10 



discontinuity line can oscillate up and down, due to the fluctuations of the piston accel- 
eration (Fig. 1). As time goes on, this oscillating curve will "travel" to the wall and come 
back to the piston, experiencing some distortions on its way, caused by the differences 
in velocities of its points (Fig. 1). When this curve comes back to the piston again, it 
may well have "turning points" where its tangent line is vertical, or even contain vertical 
segments of positive length. This produces unwanted singularities or even discontinuities 
of the piston velocity and acceleration. The same phenomena can also occur when a 
discontinuity line of the initial density 7ro(?/,w) has a negative slope. 

3 Sketch of the argument 

Our proof of Theorem |r]T| is based on large deviation estimates for the Poisson random 
variable: 

Lemma 3.1 ( ||CLS|| ) Let X be a Poisson random variable with parameter A > 0. For 
any b > there is a c> such that for all < B < bVX we have 

P(|X- A| > B\^) < 2e~^^' 



This shows that the probabilities of large deviations rapidly decay, as they do for the 
Gaussian distribution. 



The principal step in our proof of Theorem p..l| is the velocity decomposition scheme 
described next. Let VL{t,uj) be the velocity of the piston at time t > for a random 
conflguration of particles u G Ql- Let At > be a small time increment. Then the law 
of elastic collision ( [L.ID implies 

k 

VL{t + At, uj) = {l- e) Vi(t, uj) + eY,{l- e)''-^ ■ v, (3.1) 

i=i 

Here k = k{t, At, u) is the number of particles colliding with the piston during the time 
interval (t, t + At), and Vj are their velocities numbered in the order in which the particles 
collide. 

We rearrange the formula (|3.1| ) as follows: 

k 

VL{t + At, cu) = (1 - ek)VL{t, uj)+eY.Vj+ x^^^ + X^'^ (3.2) 

i=i 

where 

X^^) = VL{t.^m-ef-l + ek] 

and 
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Let us assume that the fluctuations of the velocity Vl{s,lj) on the interval (t,t + At), 
are bounded by some quantity 6V: 

sup \VLis,uj)-VLit,uj)\<SV (3.3) 
se{t,t+At) 

Consider two regions on the x, v plane: 

f/ X V -VL(t,iu) - (sgnv)5V 1 , , 1 , , 

[ x — XL{t,uj) At J 

and 

[/ X V -VL(t,Lu) + (sgnv)6V 1 , , 1 , , 

D2 = < [X, V) : — — r < ^min < < ^max } (3.5) 

[ x — XL[t,uJ) At J 

Each of them is a union of two trapezoids Di = U , i = 1,2, where denotes 
the upper and Df the lower trapezoid, see Fig. 2. 

Note that Di C -D2- The bound (|3.3|) implies that all the particles in the region Di 
necessarily collide with the piston during the time interval {t,t + At). Moreover, the 
trajectory of every point (x, v) G Di hits the piston within time At. The bound (|3.3| ) 
also implies that all the particles actually colliding with the piston during the interval 
{t, t + At) are contained in D2. 

Let us denote by the number of particles in the regions Df" for r = 1, 2 at time t. 
We also denote by k~ the number of particles actually colliding with the piston "on the 
left" , and by that number "on the right" (of course, k~ + k~^ = k). Due to the above 
observations, A;f < fc^ < fc^. 

Now, suppose that t + At < tiL. Then we show that for typical configurations u) the 
particles in each domain Dr, r = 1,2, have never collided with the piston before. There- 
fore, their number, k^, r = 1,2, satisfies the laws of Poisson distribution, in particular, 
the large deviation estimate in Lemma 3A. applies. This gives the bound (for typical uj) 

- Akf < fc± < + Ak^ (3.6) 

where 

\f = E{kf) = j pl{x,v) dx dv 



and Fjf * corresponds to the action of F~'^ = F~^/^ in the original time-space coordinates 
X, t. The deviations Akf in ( p. 6] ) can be adjusted by using Lemma |XT|. The difference 
— Af is estimated by 

A^ - A^^ = / pl{x, v) dx dv < const ■ 5V At 

Jf-\d^\d^) 

By putting all these estimates together we get tight bounds on k in (^]2|). Similarly we 
get bounds on Yl!j=i vj in (|3.2|). The following is the final result of this analysis: 

Vi(t + At, iu) - VLit, to) = V{t, uj) At + X3 (3.7) 
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Here 

V{t, u) = a [QoVlit, u) - 2Q,VL{t, lu) + Q^] (3.8) 

and Qo, Qi, Q2 are defined similarly to (|2.15|) - (|2.17|) , in which Y{t) must be replaced by 
the actual piston position XLit^u) / L. The error term in (|3.71 ) is bounded by 



, , InLVAt , , 

IXal < const (3.9) 

which corresponds to Brownian motion-type random fluctuations. 

The term T>{t,uj) in (|3.7|) represents the main ("deterministic") force acting on the 
piston. The term xs describes random fluctuations of that force. When the piston 
velocity stabilizes, then the main force V should vanish, and an "equilibrium" velocity 
ViitjUj) will be established. The latter is the root of the equation V{t,uj) = 0, which is 



t7 f+ \ ^^i - V <5i - Q0Q2 

VL{t,u}) = (3.10) 

vo 

The reason why Vl(^, ^^) converges to Vl is that T> is almost proportional to Vl^Vl, i.e. 
< i^i < V/{Vl — Vl) < E2 < 00 for some constants Ei,E2. In fact, Vf,(t, cu) is a very 
slowly changing function of t, whose derivative is small: \dVL{t,uj)/dt\ < const ■ L^^Eq. 
As a result, Vl will always stay close to Vl, more precisely 

\VL{t,uj) - VL{t,uj)\ < const ■ L"MnL (3.11) 

on the entire zero-recollision interval < t < tiL. 

Now, the piston coordinate Yl{t,uj) = Xl{tL,uj)/ L is the solution of the differential 
equation 

Yl = Vl = Vl + X4 



where \xi\ < const ■ L ^InL by ( ^.11| ). On the other hand, the deterministic piston 



coordinate Y(t) is the solution of the equation Y = W, and both V and W are given 
by the same radical expression, cf. ( p^.20D and ( |3.1(JD . Lastly, a simple application of 
Gronwall's inequality completes the proof of Theorem |1.1| on the zero-recollision interval 
(0,ri). 

The proof on the one-recollision interval (ti, t^) goes along the same lines. One major 
difference is that the number of particles k^, r = 1, 2, in the domain constructed in 



the velocity decomposition scheme is no longer a Poisson variable, so Lemma |3.1| does 
not apply directly. 

To handle this new situation, we pull the domain Df" back in time, as we did before. 
But now that puUback involves one interaction with the piston (corresponding to the flrst 
collision of the particles in with the piston, which occurs during the zero-recollision 
interval < t < tiL) . 

Since the piston position and velocity at the moment of that flrst collision are random, 
the preimage of Df" will be a random domain. Its shape will depend on the piston velocity 
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V{t,uj) during the zero-recollision interval < t < tiL. We observe that the boundary 
of the preimage of is described by a random, yet Holder continuous function, and its 
Holder exponent is 0.5 due to ( |3.9| ). Then we pick a small d > and construct a (i-dense 
set in the space of all Holder continuous functions in the spirit of a work by Kolmogorov 
and Tihomirov |[K'1| . The elements of that d-dense set can be used to construct a 
finite collection of (nonrandom) domains, so that one of them will approximate the 
(random) preimage of our Df" (we need to select the small d > carefully to ensure 
sufficient accuracy of the approximation). Now the number of particles in our random 
domain (the preimage of D^) can be approximated by the number of particles in the 
corresponding nonrandom domain. The latter has Poisson distribution, and finally we 
can apply Lemma |3]l|. This trick gives necessary estimates on k^. 

A full proof of Theorem |L1] is given in ||CLS|| . At present, we do not know if this 
theorem can be extended beyond the critical time r^,, this is an open question. Some 
other open problems are discussed in the next section. 



4 Discussion and open problems 



1. The main goal of this work is to prove that under suitable initial conditions random 
fluctuations in the motion of a massive piston are small and vanish in the thermodynamic 
limit. We are, however, able to control those fluctuations effectively only as long as the 
surrounding gas particles can be described by a Poisson process, i.e. during the zero- 
recollision interval < r < ri. In that case the random fluctuations are bounded by 
const-L~^ InL, see Remark 2 after Theorem |1 . 1| . Up to the logarithmic factor, this 



bound is optimal, see |CE] and earlier estimates by HoUey |0, Diirr et al. 

During the one-recoUision interval ti < t < T2, the situation is different. The proba- 
bility distribution of gas particles that have experienced one collision with the piston is 
no longer a Poisson process, it has intricate correlations. We are only able to show that 
random fluctuations remain bounded by see again Remark 2. Perhaps, our bound 

is far from optimal, but our numerical experiments reported in [|UL|| show that random 
fluctuations indeed grow during the one-recollision interval. 

We have tested numerically whether random fluctuations remained small after more 
than one recollision, i.e. at times r > T2. We found that for some initial ttq they actually 
increased very rapidly, and we conjectured that the rate of increase was exponential in 
T. We found, indeed, that at times r ~ logL the fluctuations became large even on a 



macroscopic scale, and then many unexpected phenomena occurred [CL 



Interestingly, the exponential growth of random fluctuations seems to be related to 
the instability of our hydrodynamical equations. We found that small perturbations of 
the initial density ttq can grow exponentially in r under certain conditions, matching the 
growth of random fluctuations of the piston motion in the mechanical model. We refer 
the reader to | |CL| ] for further discussion and to our work in progress (JCLP|. 

2. It is clear that in our model recollisions of gas particles with the piston have a very 
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"destructive" effect on the dynamics in the system. However, we need to distinguish 
between two types of recolhsions. 

We say that a recolhsion of a gas particle with the piston is long if the particle hits a 
wall X = or X = L between the two consecutive collisions with the piston. Otherwise 
a recolhsion is said to be short. Long recolhsions require some time, as the particle has 
to travel all the way to a wall, bounce off it, and then travel back to the piston before it 
hits it again. Short recolhsions can occur in rapid succession. 

We have imposed the velocity cut-off (P4) in order to avoid any recolhsions for at least 
some initial period of time (which we call the zero- recolhsion interval). More precisely, 
the upper bound Wmax guarantees the absence of long recolhsions. Without it, we would 
have to deal with arbitrarily fast particles that dash between the piston and the wall very 
many times in any interval (0, r). On the other hand, the lower bound v^i^ was assumed 
to exclude short recolhsions. 

There are good reasons to believe, though, that short recolhsions may not be so de- 
structive for the piston dynamics. Indeed, let a particle experience two or more collisions 
with the piston in rapid succession (i.e. without hitting a wall in between). This can 
occur in two cases: (i) the particle's velocity is very close to that of the piston, or (ii) 
the piston's velocity changes very rapidly. The latter should be very unlikely, since the 
deterministic acceleration of the piston is very small, cf. Theorem |2.1| c. In case (i), the 
recolhsions should have very little effect on the velocity of the piston according to the 
rule (|1.1|), so that they may be safely ignored, as it was done already in earher studies 

m,mu. 



We therefore expect that our results can be extended to velocity distributions without 
a cut-off from zero, i.e. ahowing Vmm = 0. 

3. In our paper, L plays a dual role: it parameterizes the mass of the piston (M ~ L^), 
and it represents the length of the container (0 < x < L). This duality comes from our 
assumption that the container is a cube. 

However, our model is essentially one- dimensional, and the mass of the piston M and 
the length of the interval < x < L can be treated as two independent parameters. In 
particular, we can assume that the container is infinitely long in the x direction (so, that 
L is infinite), but the mass of the piston is still finite and given by M ~ L^. In this case 
there are no recolhsions with the piston, as long as its velocity remains small. Hence, our 



zero-recollision interval is effectively infinite. As a result. Theorem |Tri| can be extended 
to arbitrarily large times. Precisely, for any T > we can prove the convergence in 
probability: 



and 



P sup \YLiT,uj) -Y{t)\ <CT\nL/L] ^ 1 

\0<T<T J 



P{ sup \WLiT,uj) -W{t)\ <CT\nL/L] ^ 1 

\0<T<T J 



as L ^ oo, where Ct > is a constant and Y{t) and W{t) = Y{r) are the solutions of 
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the hydrodynamical equations described in Section 2. 

4. Along the same hnes as above, we can assume that the container is (i-dimensional with 
d > 3. Then the mass of the piston and the density of the particles are proportional to 
L'^^^ rather than L^. 

When d is large, the gas particles are very dense on the x, v plane. This leads to a 
much better control over fluctuations of the particle distribution and the piston trajectory. 



As a result. Theorem p..l| can be extended to the fc-recollision interval {Tk,Tk+i), where 
k > 1 depends on d. It can be shown that for any k > 1 there is a c/^ > 3 such that 
for all d > dk the convergence (|1.1CI|) and (|1 . 1 1|) holds with = Tk- Therefore, a higher 



dimensional piston is more stable than a lower dimensional one. 

It would be interesting to investigate other modifications of our model that lead to 
more stable regimes. For example, let the initial density vro(?/, v) of the gas depend on the 
factor a = eL"^ in such a way that 7io{y, v) = a~^p{y, v), where p{y, v) is a fixed function. 
Then the particle density grows as a — * 0. This is another way to increase the density 
of the particles, but without changing the dimension. One may expect a better control 
over random fluctuations in this case, too. 
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Figure 1: A horizontal discontinuity line (bottom) comes off the piston as an oscillating 
curve (top). 





Figure 2: Region Di is bounded by solid lines. Region D2 is bounded by a dashed line. 
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